Localization of a Bose-Fermi mixture in a bichromatic optical lattice 
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We study the localization of a cigar-shaped super-fluid Bose-Fermi mixture in a quasi-periodic 
bichromatic optical lattice (OL) for inter-species attraction and intra-species repulsion. The mixture 
is described by the Gross-Pitaevskii equation for the bosons, coupled to a hydrodynamic mean-field 
equation for fermions at unitarity. We confirm the existence of the symbiotic localized states in 
the Bose-Fermi mixture and Anderson localization of the Bose component in the interacting Bose- 
Fermi mixture on a bichromatic OL. The phase diagram in boson and fermion numbers showing the 
regions of the symbiotic and Anderson localization of the Bose component is presented. Finally, the 
stability of symbiotic and Anderson localized states is established under small perturbations. 
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I. INTRODUCTION 



Anderson localization of super-fluid atomic gases in 
weak disordered potentials with a large exponential tail is 
currently attracting a lot of interest in both experimental 
and theoretical studies. In experimental studies, disorder 
laser speckle [l| and quasi-periodic optical lattice (OL) 
[2| have been used to localize Bose-Einstein condensates 
(BEC). The original description of Anderson localization 
was based on the localization of non-interacting quan- 
tum waves in disorder potentials due to a cancellation of 
wave fronts coming from different locations of the disor- 
der potential. However, super-fluid atomic gases usually 
are interacting and the study of Anderson localization 
has also been extended to the case of localization of a 
single-component 13 or binary Q BEC under repul- 
sive effective interactions. Both quasi-periodic OL in one 
d, Q and three dimensions and random potential Q 
were used in these studies. The effect of temperature on 
Anderson localization has also been investigated @. 

Another topic of current interest is the problem of 
trapped binary super-fluids, where the two components 
could be two different hyperfine states of the same atom 
or two different atoms. In the latter case, 
degenerate Bose-Fermi *^Rb-'"'K mixture [ll|, degener- 
ate Fermi- Fermi mixture of 12| and of ^Li were 
considered among others. In particular, the Feshbach 



resonance technique driven by a magnetic |14l | or optical 
[Tsl field allows one to vary the atomic interaction that 
opens the possibility for the study of localization of bi- 
nary condensates with controllable interactions. A Bose- 
Fermi mixture can exhibit quite distinct features as the 
number of atoms and the inter-species and intra-species 
interactions arc varied. For example, the intra-species 
repulsive interaction combined with the inter-species at- 
traction may give rise to the symbiotic soliton [l6l . [l7| . 
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In reference [l8|, we demonstrated that the inter-species 
attraction contributes to an attractive point-like effective 
potential which affects the atom density distribution and 
the stability properties of the condensate [l^ • 

Here we investigate the symbiotic localization of a 
cigar-shaped Bose-Fermi superfluid mixture without a 
trap. The Fermi superfluid is taken to be a mixture 
of equal number of spin up and down components in 
a paired state. We also consider Anderson localization 
of the Bose component in the Bose-Fermi mixture in a 
quasi-periodic bichromatic OL when the Bose component 
is not localized in the absence of the OL. In both cases 
the intra-species Bose interaction is taken to be repulsive 
and the Fermi component is considered to be at unitar- 
ity [20| . The unitarity limit is the limit of strong at- 
traction between spin-up and -down fermions with the 
Fermi-Fermi scattering length approaching infinity. Con- 
sequently, the properties of the Fermi superfluid in this 
limit become universal and independent of the Fermi- 
Fcrmi scattering length and are solely determined by the 
Fermi energy and Fermi momentum. The unitary limit of 
Fermi superfluid can now be routinely achived and stud- 
ied in laboratory [2^ . For a theoretical description of the 
interacting Bose-Fermi mixture we consider the mean- 
flcld Gross-Pitaevskii (CP) equation of bosons coupled 
to a mean-field hydrodynamic density-functional equa- 
tion (2]| for fermions at unitarity. The present investiga- 
tion involves numerical simulation as well as analytical 
study based on a variational approach. The formation 
of symbiotic and Anderson localizations as a function of 
Bose and Fermi number is illustrated in a phase diagram. 
By numerical simulation of the coupled mean-field model 
using the split-step Fourier spectral method, we confirm 
the existence of the symbiotic localized states in the Bose- 
Fermi mixture without external trapping potential. We 
also identify Anderson localization of the Bose compo- 
nent with an exponential tail assisted by the Fermi com- 
ponent in the presence of a quasi-periodic bichromatic 
OL. We find that the Anderson localization of the Bose 
component is possible in the mixture for an inter-species 
attraction which neutralizes mostly the bosonic repul- 
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sion. Finally, we establish the stability of both symbiotic 
and Anderson localizations under small perturbations. 

In Sec. ini we present a brief account of the coupled 
mean-field model and the bichromatic OL potential used 
in the study. The analytical expressions for the width 
of the Bosc and Fermi localized states obtained by the 
variational analysis of the mean-field model are in reason- 
able agreement of the numerical solution. The numerical 
results for the static and stability properties of the sym- 
biotic and Anderson states are presented in Sec. IIIII In 
Sec. IIVI we present a brief summary. 



II. ANALYTICAL CONSIDERATION OF 
LOCALIZATION 

We consider a binary super-fluid Bose-Fermi mixture 
of A'b bosons of mass me and A^f fermions of mass mp 
at zero temperature. The spin-half fermions are taken 
to be at unitarity and populated equally in spin-up and 
down states. The intra-species Bose interaction is taken 
to be repulsive and the inter-species Bose-Fermi interac- 
tion is taken to be attractive. Experimentally, this situ- 
ation is, for example, accessible in a Bose-Fermi mixture 
[13, mi of bosonic ®^Rb atoms in the hyperfine state \F = 
2, m-p = 2) and fermionic ^*^K atoms in the two equally 
populated hyperfine states \F = 9/2, mp = —9/2) and 
\F ~ 9/2, TOp = —7/2). Theoretically, the bosons arc 
treated by the usual mean-field GP equation with con- 
tact interaction, which is equivalent to the standard hy- 
drodynamical equations. The fermions are treated by a 
mean-field density-functional equation, which is equiva- 
lent to the standard hydrodynamical equations for the 
Fermi superfluid. The bosons and fermions are assumed 
to interact by a zero-range potential. Within the frame- 
work of the density-functional theory, a coupled GP-type 
equation for the mixture can be obtained pll . [23j . 

The system is made effectively one-dimensional (ID), 
assuming that the mixture is confined in transverse direc- 
tions by a tight axisymmctric harmonic potential of fre- 
quencies Lij±B and loa^f for bosons and fermions, respec- 
tively. The three-dimensional equations of Ref. [21| for 
a cigar-shaped superfluid Bose-Fermi mixture can then 
be reduced to an effective ID form by integrating out 
the dependence on the radial spatial variables [23|. The 



dynamics of the ID mixture is described by the dimen- 
sionless coupled time-dependent nonlinear equations [23j 



.duB 

i du-p 
2~dr 



gB\UB\ UB 



i d^uB 
"2 dx^ 

-gBpN^lupl'^UB + V{x)uB, (1) 
-gBFNB\uB\'^UF + V{x)uf, (2) 



/^^ \uB\^dx = /^^ luppdx = 1. In Eqs. © and ©, 
we set the transverse oscillator lengths a±B = cl±f = a_L 
and uj^B = (^±F = measuring length, energy and 
time in units of a±, huj± and i^J^, respectively [23j . 
This implies that 2mF = ruB, a condition which is 
roughly satisfied by the *^Rb-^°K mixture. The dimen- 
sionless interactions are [2^ gB = 2(aB/a±)-/VB, 5bf = 
6iaBF/a±),gF = (37r2)2/3(3^/5)iv2/3^ ^here or is the 
Bose scattering length, gbf is the Bose-Fermi scattering 
length, the universal Bertsch constant ^ has the value 
^ — 0.4 at unitarity [13, [111- We consider the scatter- 
ing length for collision between ^''Rb atoms to be posi- 
tive, Ob ~ lOSoo with oq the Bohr radius [1^, whereas 
the inter-species scattering length between bosons (f^Rb) 
and fermions (^°K) is negative cbp ~ — 284ao [H, [l^l- If 
the Bosons are decoupled from the Fermions (by consid- 
ering qbf = 0), they satisfy the GP equation and it was 
found in Ref. [j| that Anderson localization is destroyed 
due to Bose-Bose repulsion while the stationary states 
are exponentially localized [131 ■ Once a Bose-Fermi at- 
traction above a critical value is introduced via a nonzero 
negative abf, so that the Bose-Bosc repulsion in the GP 
equation is nearly compensated for by interatomic attrac- 
tion, localization is restored. 

The coupled equations (P) and will be used to study 
a cigar-shaped, localized Bose-Fermi system [28l - l30j . We 
study the symbiotic localized states of this system and 
Anderson localization of bosons in the Bose-Fermi mix- 
ture. To study the symbiotic localized states due to the 
strong inter-species attraction we take V{x) = 0. To 
study the Anderson localization, following the experi- 
ment of Roati et al. [2|, the potential V{x) is taken to 
be a quasi-periodic bichromatic OL of incommensurate 
wave lengths: 



V{x) ^^Aisiv?(kix), 
1=1 



(3) 



with Al ~ 2'K'^si/Xf, {I = 1,2), where A;'s are the wave- 
lengths of the OL potentials, s/'s are their intensities, and 
ki = 2'k/Xi the corresponding wave numbers. Without 
losing generality, we take si = 10, S2 = 0.3si, Ai — 10 
and fcs/fci = (v^- l)/2 [3l| which roughly represent the 
generic experimental situation [^ . 

To obtain the stationary localized states of the cou- 
pled equations (HJ and ([2]), we may set UB,Fix,t) = 
(/)B,p(a;) exp(— j/XB,pi) with fiB.F the respective chemical 
potentials. The real wave functions, (/)B,p(a;), obey the 
stationary equation, 



/^p0p/2 



-(I)b/2 + gB(t>l + gBFNF(t>l(t>B + 1^0B, (4) 



- gF(f>l/^ 



where V{x) is the trapping potential acting on both the 
Bose and Fermi components, ub = UBix,t) and Mp = 
Mp(a;,t) are the ID wave functions with normalization 



where the prime denotes space derivative. The stationary 
localized states can be investigated by the Gaussian vari- 
ational approach (32j . In this approach, the Lagrangian 
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density for Eqs. (|4]) and ([5]) is 

^ - 7VB[MB0|-(0i5)V2-5B4/2-V'0|] 

-5BFiVBiVF0|0^. (6) 

We use the variational trial function 



ri/4' 



wb,f 



■ cxp 



2wl 



(7) 



where the parameters i(Jb,f are the widths of the Bose 
and Fermi localized states, respectively, with normaliza- 
tion A/b.f = f(^)'^^- The trial function ([7]) and 
potential ([3]) lead to the the effective Lagrangian 
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^B.F = — 



(8) 
(9) 



Euler-Lagrangc equations, dL/dfiB.F = yield: A/b = 
A/p = 1- The remaining Euler-Lagrange equations 
dL/dwB.F = yield, respectively, 

1 + ^B^ _ 2Wb + 2ArFGB = 0, (10) 



/2tt 

TC /o \ 3/2 4/3 

1 + ^-8W^F + 8iVBGF = 0, (11) 



Wb,f = WB^F^^ikfe 



2-kfw 



2 2 



G 



5bfu'b,f 



B,F 



/7r(w| +w|) 



2x3/2' 



(12) 
(13) 



where we have set Nb = Af = 1- Equations ([TU)) and 
([TT|) give the spatial widths of the localized states. 

To study how widths u'b,f evolve with the particle 
numbers A^b.f we solve Eqs. (fTO|) and (fTT|) and show the 
phase diagram of the number of bosons and fermions in 
Fig. [1] illustrating the regions where the localization can 
exist with or without trapping potential (j3|). The region I 
bounded by the solid line corresponds to V{x) = where 
Eqs. ([TU)) and pT|) have finite real solutions for both 
widths wb.f- In this region, symbiotic localized states 
supported by inter-species Bose-Fermi attraction can be 
produced in spite of intra-species Bose and Fermi repul- 
sion. Generally, a self-repulsive Bose or Fermi super-fluid 



FIG. 1: (Color online) The phase diagram of Bose and Fermi 
atom numbers showing the regions of localization without ex- 
ternal trap (region I) and with bichromatic OL @ (regions I 
and II) from a solution of the variational equations (|10|) and 
(lllf) . In regions III and IV the variational equations have no 
solution. Single-peaked localization in the bichromatic OL is 
possible from a numerical solution of Eqs. ([1} and Q in the 
region bounded by the dashed lines. 



cannot support a localized state by itself, however, a suf- 
ficiently strong inter-species attraction can induce a net 
effective attraction responsible for the formation of sym- 
biotic states in Bose-Fermi system [T6j . 

With bichromatic OL (O, Eqs. ^ and ^ yield 
real solutions for both widths wb,f in regions I and II 
bounded by the lines with crosses. In regions III and 
IV, however, finite real solutions of Eqs. ([TU)) and ([TT)) 
do not exist. In region II, the fact that the widths wb,f 
are finite and real implies that the localized states are 
created by the quasi-periodic bichromatic OL ([3]). The 
parameter space in Fig. [T] where localization can exist is 
enlarged from region I to region II because of the bichro- 
matic OL. The phase diagram also shows that the effect 
of the bichromatic OL is small when A^b is small and 
A^F > 200. For example, the region II between regions I 
and III is very narrow. This implies that, compared to 
the bichromatic OL, the inter-species interaction is dom- 
inant in this region, because a larger A^f and small Nb 
induce a large enough effective inter-species attraction 
to localize the Bose-Fermi mixture. When A^b is larger, 
however, the effect of the bichromatic OL is more im- 
portant on localization as indicated by a wide region II 
between regions I and IV. 

The present discussion based on a Gaussian ansatz for 
the wave function has its limitation in the presence of 
the bichromatic OL where it is possible to have multi- 
peak density for the components not taken care of by the 
simple Gaussian ansatz. Single-peaked Gaussian-type lo- 
calization is possible in regions I and II. However, in the 
presence of bichromatic OL, localized states can exist be- 
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yond regions I and II into regions III and IV, where at 
least one of the components has a multi-peak density dis- 
tribution along the OL. Such states are not obtainable 
from a Gaussian variational analysis and will be studied 
in Sec. |lll] using the full numerical solution of Eqs. ([T]) 
and 




III. NUMERICAL RESULTS 



We perform the numerical integration of coupled GPEs 
dl]) and ([2]) employing real- and imaginary-time propaga- 
tion using the split-step Fourier spectral method with 
space step 0.04, time step 0.001. The time evolution is 
continued till convergence. We also checked the accuracy 
of the results by varying the space and time steps and 
the total number of space and time steps. Although the 
imaginary-time propagation method could find some of 
the localized states, the stability of these states were con- 
firmed through the real-time propagation method. For 
studying the dynamics we used real-time propagation 
corresponding to the solution of the full time-dependent 
Eqs. © and da]). 

We first confirm the existence of the symbiotic local- 
ized states. By numerical integration of GPEs ([1]) and 
([2]) with V{x) = 0, we find that the symbiotic localized 
states of the BE systems can exist in Region I of Fig. 
[TJ In order to investigate the effects of the number of 
atoms of the two components on the localized states, the 
typical numerical (N) and variational (V) widths of the 
atom density profiles of the stationary symbiotic local- 
ized states are exhibited in Fig. [5] The numerical widths 



are calculated via 



^B.F 



+ 00 



We find 



that the numerical results arc in good agreement with 
the variational results in the central part of region I. Near 
the edges (the solid line in Fig. [ij , the difference between 
the numerical and variational widths is larger because of 
the deformation of the wave functions, while the atom 
density distributions deviate from the single-peak Gaus- 
sian shape and assume a multi-peak structure. When 
A'b.f arc smaller, wc find that the deformation of the 
fcrmionic component is larger than that of the bosonic 
component consistent with Fig. [l](a). On the contrary, 
the deformation of the bosonic component is larger than 
that of the fermionic component when 7Vb,f are larger 
consistent with Fig. [5] (b) . 

A careful analysis of the density distribution of the 
localized Bose-Fermi system reveals some interesting fea- 
tures. In fact, an attractive inter-species interaction may 
induce a spatially localized effective attraction [l8| which 
can reduce the width of the atom density profiles, and 
even develop a sharp peak on the top of the bosonic and 
fermionic atom density envelopes [l9[ as we shall see in 
the following, viz. Fig. [31 Once a peak appears in the 
density of the first species (Bose or Fermi) , this induces a 
highly localized effective interaction in the second species 
(Fermi or Bose) resulting in a highly localized structure 
in the second species. The small width of a species, say 



FIG. 2: (Color online) Dimensionless numerical (N) and vari- 
ational (V) widths wb,f of the Bose and Fermi components 
of the stationary symbiotic localized states with V{x) — vs. 
(a) iVp for A^B = 1000, and vs. (b) Nb for Np = 280. 



Bose, may, however, emerge in two ways: (i) a large over- 
all attractive inter-species interaction gBpA^FlwppUB due 
to a large number of the second species, Np, or (ii) a 
highly localized inter-species interaction gBpA^FlwFpWB 
due to a narrow localized state of the second species, 
Fermi. 

Figure [2] (a) indicates that wb,f decreases monotoni- 
cally as Np increases. This is because a larger A^p induces 
a lager effective attraction among bosons [l6[ , which re- 
duces the Bose width wb due to possibility (i) above. 
Simultaneously, the narrower Bose localization causes a 
stronger attractive point-like effective potential for the 
Fermi component, which reduces the Fermi width w-p 
due to possibility (ii). Figure [3] (b) shows that, as A^b 
increases, wp decreases and wb increases. With the in- 
crease of A^B for a fixed A^p, the Bose species becomes 
more repulsive due to an increase of intra-species repul- 
sion for a roughly fixed inter-species attraction thus in- 
creasing the width of the Bose species. However, the 
increase of A^b for a fixed A^p increases the intcr-specics 
attraction on the Fermi species thus reducing the Fermi 
width. 

Next, we study the Bose- Fermi localized states with 
quasi-periodic OL ^ in region II of Fig. [TJ Typical 
atom density profiles of stationary Bose and Fermi lo- 
cahzed states are shown in Figs. [3| (a) — (d). Wc find 
that, because of the bichromatic OL, the localized Bose 
or Fermi states may exist beyond region I. For example, 
Fig- El (a) shows both Bose and Fermi components in 
single-peaked configuration for parameters corresponding 
to region II. The numerical simulation indicate, however, 
that both components can be single-peaked only in the 
region between the dashed lines in Fig. [1] Wc name this 
region the single-fragment region. In the region below 
the single- fragment region, (viz.. the region with smaller 
Ab,) the atom density profile of the Bose component is 
single peaked, but the Fermi component shows a sym- 
metrical multi-peak structure with a pronounced peak at 
center, as shown in Fig. [3] (b). However, in the region 
above the single- fragment region, (viz., the region with 
larger A^b,) the density profile of the Fermi component is 
single-peaked and that of the Bose component is multi- 
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FIG. 4: (Color online) Numerical density \ub\^ of the Bose 
component in the Bose-Fermi mixture on bichromatic OL Q 
vs. dimensionless x in logarithmic scale for (a) A'f = 90 and 
Nb = 600 (solid line), 700 (dashed line), 800 (dotted line), 
and for (b) Nb = 600 and A''f = 70 (solid line), 80 (dashed 
line), 90 (dotted line). The crosses are exponential fit to den- 
sity tails ~ exp( — |x|/Lioc), where Lioc is the dimensionless 
localization length. 
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FIG. 3: (Color online) Numerical densities |ub,f|^ of the cou- 
pled Bose-Fermi localized states with the bichromatic OL Q 
vs. dimensionless position x for (a) A''b = 800, Nf = 80, 
(b) A^B = 200, iVp = 50, (c) Nb = 1800, iVp = 90, (d) 
A^B = 1800, A^F = 90 in logarithmic scale. The solid line 
corresponds to bosons and the dashed line to fermions. 



peaked as shown in Figs, ^c) and (d). To show this 
clearly. Fig. O^d) is plotted in a logarithmic scale which 
shows that the Bose density profile develops undulating 
tails occupying many OL sites. When A^b is larger, as 
in Figs. [3] (c) and (d), the Bose system becomes more 
repulsive and hence spreads to many OL sites and thus 
develops undulating tails. However, a large A^b also in- 
troduces a strong inter-species attraction on the Fermi 
component, thus confining it in a small region. This nar- 
row Fermi peak creates a smaller narrow peak in the Bose 
component due to inter-species attraction as seen in Fig. 
[3] (c) . The role of Bose and Fermi is interchanged for a 
small Nb as can be seen in Fig. [3] (b), where one has 
a single-peak Bose distribution on top of a multi-peak 
Fermi distribution. The situation for an intermediate 
value of Nb is shown in Fig. [3] (a) . 

Now we consider the possibility of Anderson localiza- 
tion of the bosons with exponential tail in density in the 
interac ting Bose-Fermi mixture on bichromatic OL ([3]) 
[IHliSlllI- There are domains around region II of phase 
diagram shown in Fig. [T] where the inter- and intra- 
species interaction on the bosons approximately "can- 
cel" each other. Consequently, the bosons behave like 
quasi-free particles (with effective weak repulsion) and 
the weak quasi-periodic bichromatic OL is necessary for 
localizing the bosons. In order to understand the novel 
phenomenon, we investigate the localization of Bose and 
Fermi components in the bichromatic OL ([3]) away from 
the domain of strong inter-species interactions, where 
there is no localization in the absence of the bichromatic 
OL potential, e.g., around region II of Fig. [T] In that 



case, we find that the wave function of the Bose compo- 
nent may possess a pronounced exponential tail demon- 
strating Anderson localization. However, we could not 
find a similar exponential tail of fermion density distri- 
bution in the whole phase diagram in Fig. [T] 

To study Anderson localization of the Bose component, 
while both Bose and Fermi components are single-peaked 
under the action of OL trap (jSj, we plot in Fig. |4]the 
atom density |ubP, in log scale, of the stationary Bose 
localized states. The central part of the atom density 
is quasi-Gaussian, however, with a long exponential tail. 
We have also shown in Fig. |4] the exponential fitting to 
density tails with ^ exp(— |a::|/Lioc), where Lioc is the 
localization length. (Note that our definition of local- 
ization length 0] differs from that of Refs. P, Q by a 
factor of 2.) The length L\oc depends on the effective 
nonlinearity of the bosonic component and the disorder 
of the quasi-periodic OL [33j . Here, we fix the quasi- 
periodic OL ^ and study the effect of nonlinearity on 
localization length. In Fig. |l](a), for N-p = 90, Ljoc is 
= 0.25 for Nb = 600, = 0.28 for Nb = 700 and = 0.35 for 
A^B = 800 and the localization length increases with Nb ■ 
This is because that a larger A^b, for a fixed Np implies a 
larger repulsive bosonic nonlinearity, thus resulting in a 
larger value of localization length. Figure |3](b) indicates 
that, for a fixed Nb = 600, Lioc decreases with Np with 
Lioc = 0.44 for A^F = 70, = 0.32 for A^p = 80 and = 0.25 
for A^F = 90. A larger A^f, for a fixed Nb, implies a 
larg er inter-species attraction on the bosonic component 
[Iq . thus reducing the localization length. 

We also investigate numerically the Anderson local- 
ization of the Bose component ([3]) for small A^^. With 
OL (3), the Bose localization length Lioc versus A^f and 
A^B are illustrated in Figs. [5] (a) and (b), respectively. 
As shown in Fig. [S] (a), because of inter-species attrac- 
tion, Lioc for a fixed Nb monotonically decreases with 
the increase of A^f- For smaller Nb{— 200), there is a 
saturation of Ljoc and Anderson localization with expo- 
nential tail continues with the increase of A^f • For larger 
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FIG. 5: (Color online) Dimensionless localization length Lu 



of Bose density |itB 



(a) Nf for Nb = 200,400, and 



vs. (b) A^B for Nf = 100, 280, in the Bose-Fermi mixture on 
bichromatic OL ([3]). 
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FIG. 7: (Color online) The dimensionless center of the Bose 
state xb vs. dimensionless time t during the location oscil- 
lation of the Bose-Fermi mixture on the bichromatic OL Q 
initiated suddenly by introducing an initial momentum Pb.f 
through the transformation mb,f — > wb.f exp(ipB,Fa;). 




FIG. 6: (Color online) Numerical densities |ub,fP vs. dimen- 
sionless X and t of the symbiotic Bose-Fermi localized states 
after introducing a momentum in one of the components by 
transformation ub,f — >■ ub,f exp(ipB,Fa:). The parameters 
are A^b = 1000, iVp = 160 and (a) pb = 0.02, pp = 0, (b) 
Pb = 0,pf = 0.1. 



second component. By introducing an initial transla- 
tion to the Bose or Fermi component, the evolution of 
the atom density envelopes is presented in Fig. |6] for 
Nb = 1000, iVp = 160, and for (a) ps = 0.02, = 0, 
and for (b) pe ~ Q,Pf ~ 0.1. We find that the two 
components remain bound together and move with the 
same constant velocity although the initial momentum 
is given to only one component. The binding is caused 
by the stronger inter-species attraction. We see in Fig. 
[5] that the symbiotic Bose-Fermi localized states remain 
unchanged after perturbation, which confirms the stabil- 
ity of these states. 



Nb{= 400), Lioc is reduced rapidly with the increase of 
Np until the Anderson localization with exponential tail 
is destroyed and one has strong localization without ex- 
ponential tail. In Fig. [5] (b), we show the variation of 
Lioc with Nb for Np = 100 and 280. From the phase plot 
in Fig. [T]wc find that in both cases Anderson localization 
is rapidly destroyed as we enter the region I of Fig. [1] of 
permanent symbiotic trapping with the increase of Nb 
and the localization length is reduced rapidly. 

Now we study numerically some dynamics of the lo- 
calized Bose-Fermi states and test their stability. First, 
we consider a Bose-Fermi state in the region I of Fig. [1] 
corresponding to permanent symbiotic trapping. To do 
this, first we create the stationary localized Bose-Fermi 
mixture under appropriate conditions. Successively, at 
t = 0, we suddenly introduce a phase exp(jpB,F2;) in 
the wave function of one of the components to initiate 
a translational motion. No momentum is given to the 



Next we consider the dynamics of an Anderson local- 
ized Bose state in the Bose-Fermi mixture on the bichro- 
matic OL ^ with Nb = 800 and Np = 90 correspond- 
ing to one of the localized states of Fig. S] (a). The 
dynamics is initiated by introducing a phase exp(zpB,Fa^) 
in the wave function of one of the components to start a 
translational motion. Because of the quasi-periodic OL 
([3]), the bound Bose-Fermi localized state oscillates to- 
gether periodically about the origin in one site of the OL 
after introducing an initial momentum into one of the 
components, as shown in Fig. [7] which shows only the 
movement of the center of the Bose component, e.g., xb 
vs. t. Different initial momenta result in oscillation with 
different amplitudes, see the solid line (pb = 0.3), dashed 
line {pp ~ 0.3) and the dotted line (pb = 0.1) in Fig. [T] 
Sustained oscillation of the Bose (as well as the bound 
Fermi component not shown here) component confirms 
the stability of the coupled Bose-Fermi state. 
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IV. SUMMARY 

Here we studied symbiotic localization in a cigar- 
shaped Bose-Fermi super-fluid with intra-species repul- 
sion and inter-species attraction. The Fermi component 
is considered at unitarity. In the presence of a quasi- 
periodic bichromatic OL, the Bose component in the 
Bose-Fermi mixture could exhibit Anderson localization 
with large exponential tail due to a near cancellation of 
the intra-species repulsion and inter-species attraction 
while no localization is possible in the absence of the 
OL. No evidence of Anderson localization with long ex- 
ponential tail was found in the Fermi component. The 
cigar-shaped Bose-Fermi mixture is described by an ef- 
fective one-dimensional GP equation for bosons coupled 
to a mean-field hydrodynamic equation for fermions at 
unitarity [23| . In this study, we use both numerical and 
variational solution of the mean-field equation for the 
mixture. We obtain a phase plot of number of bosonic 
and fermionic atoms showing the domain of symbiotic lo- 
calization without external trap and localization in the 



presence of a quasi-periodic bichromatic OL. The numer- 
ical and variational densities of the localized states are 
in good agreement with each other. Both symbiotic lo- 
calization and Anderson localization are found to be dy- 
namically stable when given a small initial velocity to 
one of the components. The Bose and Fermi components 
are found to move together under such a perturbation. 
With present know how, it should be possible to study 
the Anderson localization of the Bose component in the 
Bose-Fermi mixture on a bichromatic OL experimentally 
under appropriate condition. 
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